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ABSTRACT 

The  damage  tolerance  assessment  of  complex  aerospace  structural 
components  requires  the  capability  of  effective  modeling  of  3D  cracks  and 
their  associated  propagation  and  velocity  and  path  under  fatigue  loads.  A 
3D  mixed  mode  crack  propagation  theory  is  presented  which  includes  the 
effect  of  Ki,  Kn,  and  Km,  as  well  as  non-proportional  loading,  elastic  and 
fracture  resistance  anisotropy,  and  fracture  mode  asymmetry  (viz.  the 
ability  to  transition  between  competing  tensile  and  shear  modes  of 
propagation).  A  modified  strain  energy  release  rate  criterion  including  the 
modeling  of  crack  closure  is  developed  and  presented  for  a  representative 
problem.  An  elementary,  mode  I  characterization  of  closure  is  used, 
leaving  shear  mode  closure  as  fertile  ground  for  further  study. 

Use  of  the  model  is  presented  for  an  example  problem  with  steady- 
vibratory  interaction. 


INTRODUCTION 


Three-dimensional  fracture  simulation  has  advanced  significantly  over  the  last  few 
decades.  Early  work  was  focused  on  building  the  framework  to  appropriately  represent 
cracks  in  complex  geometries,  and  calculate  sufficiently  accurate  mixed-mode  stress 
intensity  factors  [1].  As  crack  propagation  capability  followed,  much  of  the  effort  went 
to  the  development  of  the  framework  necessary  to  model  the  extending  crack  with 
minimal  user  workload,  such  as  the  FRANC3D  code  developed  at  Cornell  University  [2]. 
Non-planar  crack  growth  algorithms  typically  utilized  two-dimensional  mode  Eli  crack 
turning  theories  that  were  well  established  a  generation  earlier. 

This  approach  has  worked  well  with  a  wide  range  of  engineering  applications.  However, 
as  the  technology  to  model  non-planar  cracks  in  complex  geometries  has  matured,  the 
problem  set  has  become  more  demanding,  requiring  propagation  criteria  that  include  such 
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things  as  HCF/LCF  interaction,  Mixed  Non-Proportional  Loading  (MNPL),  fracture 
mode  asymmetry,  and  both  elastic  and  fracture  resistance  anisotropy. 

Legacy  crack  propagation/turning  criteria,  such  as  the  Maximum  Tangential  Stress 
(MTS)  criterion  [3],  assume  proportional  loading  (Kn/Ki  =  constant),  and  predict  crack 
growth  along  Kn  ~  0  path.  For  non-proportional  loading,  the  relative  proportions  of  Ki, 
Kn,  and  Km  vary  with  time  throughout  the  cycle,  and  there  is  no  path  that  enforces  Kn=0 
during  the  entire  cycle  (not  to  mention  the  influence  of  mode  III).  Though  widely 
neglected,  MNPL  can  result  from  any  structural  situation  wherein  steady  and  cyclic 
stresses  are  misaligned,  as  in  the  vibrating  blade  problem  illustrated  in  Figure  1 
Vibrating  Rotor  Blade  Resulting  in  MNPL  . 


Steady  Spin  Load 


ff 


Vibratory 

Load 


Figure  1  Vibrating  Rotor  Blade  Resulting  in  MNPL 


In  order  to  enable  crack  growth  simulation  for  this  class  of  problems,  two  significant 
advances  were  required. 

1.  A  validated  theoretical  approach  to  handle  crack  growth  and  trajectory  under 
MNPL  loading. 

2.  A  reduced-order  nonlinear  dynamics  approach  to  enable  FEM  vibration  modeling 
including  contact  between  opposing  crack  faces. 

The  current  paper  will  deal  exclusively  with  the  first  of  these  requirements,  though 
promising  approaches  to  the  second  problem  are  a  separate  subject  of  investigation  [4]. 
Theoretical  enhancements  have  been  implemented  in  the  FRANC3D  code  as  described 
herein. 


Theoretical  Background  for  Crack  Growth  with  Mixed  Non-Proportional  Loading 
Conditions 


Maximum  Stress  Criteria 

Resolving  the  (isotropic)  lead  crack  stress  intensities,  (Ki,  Kn,  Km),  into  the  asymptotic 
stress  intensities  (ki,  kn,  km)  associated  with  an  infinitesimal  crack  branch  at  angle  AO,  we 
obtain 
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(1) 


kj  (A 0)  =  <Jggyl2 7rr  =  cos  Kj  cos2  ^ ^  Kn  sin  AO 
kn  (A 0)  =  arg  Vl 7rr  =  cos  ^ -  [ Kj  sin  A 0  +  Kn  (3 cos  AO  - 1)] 

^  Q 

k,„  (A 0)  =  erft  =  Kin  cos  — 


For  the  coordinate  system  illustrated  in  Figure  2 


Figure  2  Reference  coordinate  system  for  theoretical  development 

A  similar  solution  for  3D  anisotropy  based  on  the  work  of  Hoenig  [5],  will  not  be  further 
discussed  here,  but  has  been  be  treated  and  implemented  in  FRANC3D  in  much  the  same 
manner  as  the  isotropic  theory  that  will  now  be  further  discussed. 

For  proportional  loading,  the  classical  Maximum  Tangential  Stress  (MTS)  theory, 
proposed  by  Erdogan  and  Sih  [2]  for  isotropic  materials,  asserts  that  the  crack  will  grow 
toward  the  location  of  the  maximum  tangential  tensile  stress  (equivalent  to  maximizing 
ki).  By  differentiating  ki  and  equating  to  zero 

Kn  _  -sin  A0C 
Kj  (3cosA6(.  - 1)  ^2) 
or  [6] 

A  6C  =  2tan~‘ 

where  A6C  is  the  kink  angle.  This  criterion  seeks  out  a  mode  I  crack  path.  As  illustrated 
in  Figure  3,  the  maximum  tangential  stress  theory  works  well  for  low  ductility  materials 
like  PMMA,  but  fails  to  predict  a  transition,  observed  for  7075-T6  and  2024-T3 
aluminum  alloys,  to  a  path  that  is  associated  with  the  Maximum  Shear  Stress  (MSS)  .  In 


l-Vl  +  SQ^/^)2 
4  (Kh/K,) 
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these  materials,  the  transition  is  sudden,  and  the  crack  seems  exclusively  dominated  by 
either  mode  I  or  mode  II. 

While  it  should  be  noted  that  the  data  in  Figure  3  is  for  cracks  loaded  quasi-statically 
until  the  crack  begins  to  tear,  as  opposed  to  cyclic  loading,  similar  behavior  can  occur  as 
a  result  of  cyclic  loading,  as  will  be  shown. 


Mode  mixity,  -ATANf K./K,) 

Figure  3  Test  data  of  various  investigators  [7,  8,  9,  10],  showing  transition  to  shear  mode 
dominated  growth  for  ductile  metals  at  high  mode  mixities 


Fracture  Mode  Interaction  and  Asymmetry 

Chao  and  Liu  [11]  describe  the  sharp  transition  behavior  as  a  result  of  competing  failure 
mechanisms  in  mode  I  and  mode  II  as  shown  in  Figure  4.  According  to  their  hypothesis, 
the  two  modes  do  not  interact  (see  also  [12]),  and  the  crack  will  fail  in  Mode  I  unless 


k,(  A<?) 


K 


1C 


< 


ki,(A0) 


K 


lie 


(4) 


Where  and  Knc  are  the  pure  mode  fracture  toughness  values  for  a  straight  growing 
crack,  and  the  subscript  “max”  denotes  maximizing  with  respect  to  the  kink  angle.  For 
the  purpose  of  further  discussion,  this  criterion  will  be  described  as  the  “Modal”  fracture 
criterion 
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Figure  4  Transition  mechanism  proposed  by  Chao  and  Liu  {11] 


Note,  however,  that  the  HY130  steel  in  Figure  3  behaves  in  a  different  manner1, 
exhibiting  a  more  gradual  transition  to  a  kink  angle  of  zero,  rather  than  the  angle 
associated  with  maximum  mode  II.  Thus,  both  tensile  and  shear  modes  of  growth  appear 
to  be  contributing  to  failure.  This  type  of  behavior  correlates  well  to  a  Modified  Strain 
Energy  Release  Rate  (MSERR)  approach  proposed  by  Kfouri  &  Brown  [13],  that 
suggests  a  failure  locus  of  the  form 

ITMY 


For  a  non-critical  load  state,  the  most  critical  crack  growth  direction  would  be  obtained 
by  maximizing  the  left  hand  side  of  this  equation  with  regard  to  6.  A  measure  of  how 
nearly  critical  the  loading  is  in  terms  of  an  equivalent  mode  I  stress  intensity  can  be 
written  by  solving  Equations  (4)  and  (5)  for  Kic  respectively  and  rewriting  so  that  the 
crack  is  critical  if  kieq  =  Kic. 

.  f.  /  .  l  Kir  l.  /.  „\l  ]  - -  ... 


=  MAx\kI(A0) 


=  JMa£)]2  + 


\k(m\ 


[U  AS)]2 


(Modal) 


(MSERR)  (7) 


1  The  HY-130  tests  exhibited  tearing  in  a  zig-zag  microscopic  shear  mode.  The  crack 
turning  angles  shown  for  this  alloy  were  read  from  Maccagno  and  Knott’s  photographs, 
and  reflect  the  average  trend  of  the  zig-zag  line,  consistent  with  our  intent.  Maccagno 
and  Knott  gave  quite  different  values  of  the  turning  angle. 
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While  the  modal  criterion  assumes  non-interaction  of  modes  I  and  II,  modes  II  and  III  are 
both  shear  modes,  and  can  be  combined  into  a  resolved  shear  stress  intensity 

^3  =  yjkj,  +  k]n  (8) 


It  seems  reasonable  to  assume  that  such  modes  would  interact  in  materials  that  fail  by 
either  the  Modal  or  MSERR  criterion,  which  infers  that  equations  (6)  and  (7)  can  be 
generalized  in  3D  to 


kIeq=MAX< 


'Ejl'2 

\Kjjc  j 


[*„Mf  + 


K 


ic 


V  K  me  J 


[*,,,  (A0)]2 


(9) 


(Modal) 

and, 


k  I  eq 


[*,MP  + 


V  K„c  J 


Ma6>)1  + 


V  K„IC  J 


[k,„  (AS)]2 


(10) 


(MSERR) 

The  ratios  Kic/Knc  and  Kic/Kmc  are  measures  of  the  fracture  mode  asymmetry  of  a 
material,  a  term  coined  by  Kfouri  describing  the  relative  fracture  resistance  of  a  given 
material  in  the  different  modes.  The  material  dependent  transition  points  observed  in  the 
data  in  Figure  3  can  be  correlated  to  different  values  of  these  ratios.  While  the  fracture 
mode  asymmetry  is  couched  in  terms  of  fracture  toughness  ratios,  the  intent  is  to  view 
these  ratios  as  material  parameters  in  their  own  right,  that  can  be  used  to  model  the 
transitions  between  tensile  and  shear  crack  propagation  under  both  monotonic  and  cyclic 
loading  conditions. 


Fracture  Resistance  Anisotropy 

In  addition  to  fracture  mode  asymmetry,  which  addresses  the  relative  fracture  resistance 
of  a  material  to  the  different  fracture  modes,  there  is  also  the  potential  for  the  fracture 
resistance  to  vary  as  a  function  of  crack  orientation  within  the  material.  Buczeck  and 
Herakovitch  [14]  expressed  the  fracture  resistance  in  2D  as  a  simple  elliptical  function. 
This  was  generalized  to  include  a  more  flexible  interpolation  function,  as  shown  in  Figure 
5,  and  extended  to  3D  [15,  16,  17]  to  express  the  (pure  mode)  stress  intensity  factor  at 
which  a  crack  will  propagate  in  an  arbitrary  orientation.  The  crack  orientation  in  3D  is 
identified  by  the  normal  and  tangential  vectors  (n,a)  as  illustrated  in  Figure  6a. 


6 


KPm 


2  ^  ^ 


cos“  0 

vW; 


+ 


2  s-i  '\ 


sin'  0 
yKp(90  °)2y 


i 

2n 


Figure  5  Polar  interpolation  function  for  fracture  resistance  anisotropy  in  2D  [16,  17] 
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/Crack 


Crack  Front 

a.  Geometry  of  crack  orientation  at  a  point  on  an  arbitrary  crack  front 


b.  Principle  fracture  resistance  values  and  traces  of  crack  growth  direction  a  in  cardinal  planes 
Figure  6  Physical  Parameters  Governing  3D  Fracture  Resistance  Anisotropy  [15] 


Assuming  three  principle  planes  of  material  symmetry,  each  with  their  own  2D  fracture 
resistance  interpolation  functions  as  shown  in  Figure  6b,  the  stress  intensity  at  which  a 
crack  will  propagate  in  the  (n,a)  orientation  is  given  by 
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Kp(  n,a)  = 


f  2  A 

n , 


+ 


z  2  V 

n, 


+ 


Z  2  A 


v^2y 


1 

21 


(ID 


Where  the  trace  fracture  resistance  values,  K,  of  the  a  vector  in  fracture  resistance  space, 
as  illustrated  in  Figure  6b,  are  given  by 


a/1-  a'i 2 

z  2  X 

aj 

m 

+ 

z  2  X 

ak 

1 

m  2  m 

W) 

k/J 

(non-repeating)  (12) 


and  the  K,,  are  the  six  principal  fracture  resistance  values  in  3D  (two  for  each  principal 
planes),  and  n,  and  a,  are  the  components  of  the  n  and  a  vectors.  Depending  on 
processing  symmetries,  the  six  principal  fracture  resistance  values  may  not  all  be  unique. 
If  they  are  all  equal,  and  £m=l,  the  properties  are  isotropic  with  regard  to  fracture 
resistance. 

For  the  purposes  of  the  current  development,  Kp  will  be  considered  to  be  normalized  to 
the  fracture  resistance  in  a  reference  orientation,  likely  chosen  to  be  the  orientation  for 
which  the  properties  are  most  fully  characterized.  By  definition,  Kp=l  in  the  reference 
orientation. 


Extension  to  Fatigue  Crack  Growth 

For  application  to  fatigue  crack  growth,  the  six  principal  fracture  resistance  parameters 
and  the  fracture  mode  asymmetry  parameters  are  assumed  to  be  the  same  for  all  crack 
growth  in  a  given  material,  whether  near  threshold  or  approaching  the  fracture  toughness. 
This  implicitly  assumes  that,  the  AKeff  vs  da/dN  curves  are  parallel  regardless  of  crack 
orientation  or  modality  of  failure,  and  any  differences  can  be  represented  by  an 
appropriate  horizontal  shift  (in  Kmax  and  AKeff).  While  this  is  almost  certainly  an 
oversimplification,  it  represents  an  attempt  to  include  all  these  real-world  effects  in  the 
simplest  way  possible2. 

When  we  speak  of  AKeff  in  this  venue,  we  are  referring  to  the  cyclic  stress  intensity  that 
actually  makes  it  to  the  crack  tip,  as  opposed  to  the  globally  applied  AK.  Among  the 
potential  crack  tip  shielding  effects,  Mode  I  plasticity  induced  crack  closure  [18],  though 
still  a  topic  of  lively  discussion,  remains  an  industry  standard  approach  to  account  for  R- 
ratio  effects.  Add  shear  modes,  and  there  is  the  potential  for  friction  to  reduce  the 
effective  shear  cyclic  stress  intensity.  While  there  is  ample  evidence  for  shear  mode 
crack  tip  shielding  in  the  literature  [19,  20,  21]  no  quantitative  theoretical  framework  has 
become  widely  accepted  for  modeling  this  behavior. 

Despite  these  difficulties,  it  was  recognized  that  without  taking  into  account  shielding 
effects  in  some  way,  known  phenomena  could  not  be  predicted.  Also,  it  was  desired  to  at 
least  maintain  industry  standard  capability  for  mode  I  problems,  including  the  ability  to 
predict  the  effect  of  R-ratio.  It  was  thus  decided  to  adopt  as  a  baseline  the  NASGRO 
crack  growth  model  formulation  [22,  23],  including  the  Newman  closure  equations  for 


2  Kfouri  suggested  that  the  fracture  mode  asymmetry  might  also  be  a  function  of  orientation,  but  that 
possibility  was  excluded  in  the  current  formulation. 
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Mode  I  closure,  and  attempt  to  extend  the  equations  in  a  rudimentary  way  to  account  for 
mixed  modes.  While  it  was  recognized  that  this  could  potentially  have  serious 
shortcomings,  it  would  serve  at  least  as  a  temporary  member  in  the  overall  theoretical 
framework,  with  the  opportunity  to  improve  upon  it  as  better  methods  become  available. 
For  the  purposes  of  the  present  formulation,  we  will  again  invoke  the  steady- vibratory 
scenario  of  Figure  1,  and  introduce  the  notation 

A K  =K  +  -K  (13) 

Superscripts  refer  to  the  extreme  values,  max  and  min,  for  Ki,  and  to  the  corresponding 
extreme  values  for  the  shear  modes.  That  is,  the  positive  sign  will  correspond  to  the 
extreme  load  state  with  the  most  positive  K|  value,  regardless  of  the  sign  or  magnitude  of 
the  shear  modes. 

Referring  to  Figure  3,  there  are  two  possible  assumptions  that  could  be  made  with  regard 
to  crack  closure: 


-Closure  occurs  in  the  lead  crack  only  (no  infinitesimal  kink,  resolved  stresses 
only) 

-Closure  behavior  occurs  in  the  infinitesimal  kink  tip 

With  the  assumption  of  kink  tip  closure,  crack  growth  is  evaluated  at  the  infinitesimal 
kink  tip,  and  Eq  13  would  written  in  lower  case  k’s  ,  and  be  evaluated  (and  the  sign 
convention  established)  at  the  kink  tip.  While  capital  (lead  crack)  notation  will  be 
followed  on  the  next  few  equations,  bear  in  mind  that  they  would  be  written  in  lower  case 
(kink  tip)  notation  for  the  second  assumption  above.  With  that  in  mind,  crack  closure 
will  be  defined  by  the  mode  I  component. 

A Kf  =  K,+  -  K°p  =  A K,Fl  (14) 


Where  Newman  and  global  closure  options  as  shown  in  Figure  7. 


Fi=< 


0 

1 

f(R) 

min(K  +  ,AK  j) 

A Kj 


if  K+  <  0 

for  Ak  applied 

for  Newman  closure 

for  global  closure 


(15) 


Figure  7  Illustration  of  Newman  and  global  closure 


The  Newman  closure  function  (slightly  modified  by  the  additional  material  parameter  C 
for  negative  R-ratios)  is  given  by 
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[s]=<=tr«l) 

L  J  A Kj  (l  -R) 


where 


Sop /S+  =Ao+CA1R  for 

Sop /S+  =  MAX(Ao+AlR  +  A2R2+A2R3,R) 


A 

A 

A 


=  (.825-. 34a +  . 05a 2 


nS^ 


2  a 


O  J 


1 

a 


-  (.415  -.071a)S+/cr, 

=  i- A  -A  -A 


=  2  A  +  A  - 1 


a  =  constraint  factor 


fl.O  plane  stress 
[3.0  plane  strain 


cr0  =  flow  stress 


<R<  0 
for  R>  0 


(16) 


While  AKu  and  AKm  do  not  “close”  in  the  sense  that  mode  I  does,  it  can  first  be 
postulated  that  once  the  kink  tip  closes  in  mode  I,  a  “stick”  or  “slip”  condition  exists, 
altering  AKn  and  AKra  as  follows 
AKf  =A  K„F1} 

A  K*  =AK,„F1, 

where  (17) 


F,  for  stick  friction 

F23  =  < 

[  1  for  slip  friction 

For  simple  sliding  friction,  F23  should  lie  between  stick  and  slip  values.  For  (Newman) 
closure  over  a  small  region  near  the  tip  of  the  crack,  the  weight  functions  for  tension  and 
shear  modes  are  identical,  thus  the  maximum  amount  of  K23  that  can  be  dissipated  in 
Coulomb  friction  by  a  compressive  Ki  is 

Kff  =  juK'fl'"1’  where  Kfmp  =  Max( 0,  Kf  -  K, )  (18) 


To  the  degree  that  the  global  weight  functions  for  tensile  and  shear  are  in  agreement, 
these  expressions  will  also  be  approximately  true  for  global  closure.  Allowing  the 
possibility  of  closure  at  both  extremes  of  the  cycle, 
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Maxi 


1- 


ju[Kfmp+  +  KcIomp- 
>  K\  +A  Kl 


for  Coulomb  friction 


(19) 


Note  that  like  the  stress  intensity  factors,  the  closure  parameters  Fj  and  F23  may  be 
written  in  lower  case  when  they  refer  to  values  evaluated  at  the  kink  tip. 
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Lead  Crack  Closure 

The  lead  crack  stress  intensity  factors  and  R-ratio  for  lead  crack  closure  may  be 
expressed  as 


1 

+1 

1 _ 

~K,~ 

~K,~ 

Kf, 

= 

K„ 

± 

K„ 

K±m_ 

_K in  _ 

St 

_Km  _ 

R,=—  (21) 

1  k; 


(20) 


The  opening  stress  intensity  factors  are  given  by 


Kf 

~Ki~ 

"  *i  (*/1 

~KJ)  " 

K°p 

''  ii 

= 

K 

- 

F23(Kh 

~Kn) 

(22) 

Kop 

/// 

K 

Ftz&m 

-Km) 

In  the  special  case  where  Ki+<0,  Ki+  and  Ki°p  must  are  set  to  zero  for  compression- 
compression  loading.  The  effective  crack  tip  Ak  values  are  calculated  from  the  lead  K+ 
and  Kop  values  using  (1). 


(23) 


Akf  ~ 

k](A6)-kop(A9) 

A  kf 

= 

abs{k;i(Ad)-k(;p(Ad)\ 

A  k«_ 

abs[kjn  (Ad)  - k°np (A#)] 

It  is  postulated  that  the  crack  will  grow  in  the  direction  corresponding  to  the  maximum 
growth  rate,  da/dN,  associated  with  these  Ak  values,  which  will  require  a  numerical 
search  to  maximize  an  appropriate  growth  parameter.  Note  that  the  use  of  resolved 
Akieff(A0)lmax  as  the  crack  turning  criterion  with  global  closure,  proportional  loading,  and 
positive  Ki  values  yields  results  equivalent  to  the  MTS  crack  turning  criterion.  A  more 
general  criterion  will  be  proposed  later  on. 

One  area  of  discomfort  evident  in  the  foregoing  formulation  is  the  treatment  of  negative 
R-ratios.  In  NASGRO  and  other  legacy  codes,  negative  K’s  are  nominally  allowed,  and 
act  to  accelerate  the  crack  growth  via  negative  R-ratios  using  the  Newman  closure 
equations  (an  approach  with  its  own  shortcomings  even  for  mode  I  loading).  In  mixed¬ 
mode  situations,  however,  allowing  negative  Ki  values  (which  unrealistically  involves 
crack  faces  passing  through  each  other)  leads  to  resolved  mode  II  components  at  non¬ 
zero  0.  Yet,  neglecting  negative  Ki  values  neglects  the  associated  acceleration  associated 
with  -R  ratios. 


Kink  Tip  Closure 

The  lead  crack  stress  intensity  factors  needed  to  evaluate  kink  tip  closure  are  the  same  as 
given  in  Equation  (20),  except  that  for  the  reasons  just  discussed,  in  the  case  of  a  negative 
R-ratio  at  the  lead  crack  it  was  found  necessary  to  enforce  global  closure  (truncate 
negative  Ki  at  zero)  on  the  Ki  term  before  evaluating  its  contribution  to  mode  II  at  the 
kink  tip.  Negative  lead  crack  Ki’s  are  allowed  to  contribute  to  ki,  to  preserve  the  negative 
R-ratio  acceleration  for  near- mode  I  scenarios. 
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Theoretically,  the  sign  convention  as  to  which  is  the  +  side  of  the  cycle  is  not  decided 
until  the  kink  tip  stress  intensity  factors  and  ranges  are  calculated  from  the  lead  crack 
values  using  (1). 


Akj 

"  kf  (Ad)  -  kf  (Ad) " 

A  ku 

=  abs 

kfj(AO)  -ku{A6) 

A  km 

km  (AO)  -  km  {AO) 

(24) 


However,  by  taking  the  absolute  value  and  by  defining  the  mode  I  R-ratio  as 


MINfkf  ,  kf ) 
MAX(kf,k 7) 


(25) 


we  rectify  the  use  of  the  lead  crack  sign  convention  (as  will  be  seen,  the  sign  of  the  shear 
mode  ranges  is  later  squared,  and  is  thus  inconsequential).  The  kink  tip  effective  stress 
intensity  ranges  then  given  by 


1 

> 

_ 1 

1 

> 

_ 1 

A  kf 

= 

f 23  Ak  ii 

1 

> 

1 _ 

f 23  Ak  in 

(26) 


Generalized  Crack  Propagation  Criteria 

As  alluded  to  earlier,  determination  of  the  crack  growth  direction  will  require  a  numerical 
search  to  maximize  an  appropriate  growth  parameter.  Following  the  approach  of  Buzcek 
and  Herakovitch  [14]  used  successfully  in  FRANC2D  [24]  for  many  years,  the  crack  is 
postulated  to  grow  in  the  direction  so  that  the  ratio  between  the  crack  driving  force  and 
the  crack  growth  resistance  is  maximized. 


Crack  driving  force{AO)  | 

[ Crack  growth  resistence{A6)\ 

Writing  Equations  (9,  10)  in  terms  of  Ake^  as  crack  driving  forces,  and  using  (11)  for  the 
crack  growth  resistance,  we  can  write  the  generalized  crack  propagation  criteria  as 


A kfeq  =  MAX 


r  A  kf  A 

v^7  , 


K„ 


K 


ic 


V  K  nc  J 


A  k 


eff " 


K 


ic 


V  K  me  J 


A  k 


eff‘ 

III 


(28) 


(Modal) 

And, 


(29) 


(MSERR) 

Note  that  these  Akieqeff  values  are  effective  in  the  sense  that  they  are  closure  adjusted,  and 
equivalent  in  the  sense  that  they  include  all  mixed  mode  effects  (to  the  extent  the  theory 
is  capable),  and  can  thus  be  used  in  conjunction  with  a  standard  mode  I  closure  model. 
The  angle  resulting  from  the  maximization,  9C,  is  the  predicted  kink  angle. 
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The  problem  is  that  Kmax  and  AKappiied,  as  required  in  Equation  (30),  are  not  defined  in  a 
manner  sufficiently  general  to  include  mixed  mode  behavior.  It  is  proposed  that  Kmax  be 
generalized  by  using  the  equivalent  value. 

*“  =MAx\k^,k;J  (34) 


Where  Aki  eq  can  be  written  for  the  two  material  behaviors  as 


(Modal) 

and, 


The  AKappiied  value  to  be  used  to  shape  the  near-threshold  regime  is  more  difficult  to 
define  with  confidence,  in  part,  because  there  is  some  evidence  [25  ]  that  as  threshold  is 
approached  (very  small  scale  yielding),  the  modal  transition  sometimes  disappears 
(suggesting  that  the  fracture  mode  asymmetry  ratios  reduce  near  threshold).  Based  on 
this  tentative  observation,  one  might  simply  assume  that  AKth  in  Equation  (30)  is 
evaluated  using  (31)  with  R=  Ri,  and  use  Aki  (calculated  with  no  closure,  thus  the  same 
as  in  (24))  as  AKappiied.  However,  this  would  require  that  9C  be  determined  by  maximizing 
(30)  instead  of  (28)  or  (29).  Such  an  approach  might  have  merit,  driving  cracks  in  a 
mode  I  direction  near  threshold,  but  could  potentially  be  non-conservative  if  shear  modes 
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were  able  to  contribute  to  propagation  in  the  near- threshold  regime.  On  the  other  hand, 
one  could  conservatively  neglect  threshold  altogether. 

For  the  current  implementation  it  was  decided  to  calculate  pure  mode  thresholds  as 
follows. 


A K 


th  Calculated  per  NASGRO 

l  model  with  R=Rj 


AiC=— AiC  (37) 

KIc 

ak*=^ak; 

KIc 

The  use  of  the  intrinsic  (closure-free)  thresholds  as  the  basis  for  estimating  shear  mode 
threshold  is  almost  certainly  conservative.  The  equivalent  applied  values  may  then  be 
calculated. 


kapp,ied  =  CR1T< 


C  applied  'N 


K 


J  AG= AG r 


K. 


r  A K'^2 
K^ftJ 


i  applied^ 
KII 


r  A Kf  ^ 


i  applied^ 
KIII 


y  AG= AG c 


(modal) 


(38) 


t  applied 
KIeq 


in 


A k 


applied ^ 


+ 


r  A K?^2 

\AkJi  j 


applied^  ^ 


A K 


th  \ 


k*K%j 


A  k 


applied ^ 


III 


AG=AGr 


(MSERR)  (39) 

AKi th  and  Aki  eqapplied  are  substituted  for  Akth  and  Akappiied  in  Equation  (30)  to  calculate 
the  crack  growth  rate.  This  method  is  not  expected  to  perform  well  at  predicting  MPNL 
threshold  behavior,  but  should  serve  as  a  lower  bound  for  the  purposes  of  crack  growth 
analyses  until  a  better  method  becomes  available. 


Correlation  with  Tension-Torsion  Data 


A  mixed  mode/non-proportional  loading  test  program  was  undertaken  to  provide  MNPL 
crack  turning  data.  Specimens  were  made  of  IN718  nickel  alloy,  machined  into  a 
tension-torsion  configuration  with  through-wall  cracks,  and  thus  were  predominantly 
loaded  in  mode  Eli.  Tests  included  specimens  with  tension  and  torsion  loading  in-phase, 
tension  constant  and  cyclic  torsion,  torsion  constant  and  cyclic  tension,  and  tension  and 
torsion  loading  180  degrees  out-of-phase.  Fabrication  and  testing  took  place  at  NASA 
MSFC,  and  fractography  and  K-solution  development  was  performed  by  Shelby 
Highsmith  [26]  at  Georgia  Tech.  A  summary  of  the  test  data  is  provided  in  Table  1. 
Some  crack  tips  were  observed  to  bifurcate  initially,  with  one  of  the  kinks  subsequently 
becoming  dominant.  In  other  cases,  a  negative  kink  would  occur  on  one  crack  tip,  and  a 
positive  kink  on  the  other  (speaking  in  crack  tip  coordinates,  so  a  “symmetric”  looking 
kink  arrangement  is  actually  one  positive  and  one  negative).  In  all  these  cases,  both 
primary  and  secondary  angles  were  recorded  (though  it  wasn’t  always  clear  which  was 
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which),  and  when  more  than  one  like  angle  resulted,  they  were  averaged  in  the  table.  For 
further  detail,  see  [26] . 


Table  1.  lnco718  Crack  Kink  angle  Data 


Measured  Kink  Angles 

Spec  No 

Rl 

Rll 

Rill 

KI+ 

KII+ 

KIII+ 

Mean  Beq 

Primary 

Secondary 

In  Phase 

1 

0.6 

0.6 

0.6 

19.160 

10.280 

5.6100 

28.22 

-37.0 

2 

0.6 

0.6 

0.6 

19.435 

10.670 

5.7000 

28.77 

-37.0 

8 

0.1 

0.1 

0.1 

10.940 

10.830 

2.1700 

44.71 

-38.0 

8.5 

9 

0.1 

0.1 

0.1 

11.025 

10.725 

2.7600 

44.21 

-56.5 

17 

12 

0.1 

0.1 

0.1 

1 6.470 

7.170 

1.7200 

23.53 

-28.0 

16 

0.1 

0.1 

0.1 

13.600 

8.850 

2.3800 

33.05 

-32.5 

Kl  Const 

3 

1 

0.6 

0.6 

17.705 

1 1 .090 

3.2300 

26.33 

-54.5 

4 

1 

-1 

-1 

17.775 

11.175 

2.5200 

0.00 

0.0 

6 

1 

0.0204 

0.0204 

20.405 

10.430 

6.3750 

13.84 

4.5 

10 

1 

0.1 

0.1 

10.445 

10.365 

2.2000 

25.22 

1.0 

13 

1 

0.1 

0.1 

16.130 

7.195 

1.0650 

13.30 

-46.5 

5.5 

Kll  Const 

5 

0.6 

1 

1 

21.470 

10.755 

7.0400 

33.23 

-16.0 

7 

0.1 

1 

1 

15.480 

15.470 

2.4250 

64.63 

-15.5 

11 

0.1 

1 

1 

16.240 

6.970 

1.9600 

50.06 

-10.5 

Out  of  Phase 

17 

0.1 

10 

10 

14.130 

0.902 

0.2770 

42.37 

-38.0 

18 

0.1 

10 

10 

9.960 

0.964 

0.2770 

44.81 

-64.0 

6.0 

19 

0.1 

10 

10 

16.075 

0.691 

0.1965 

39.68 

14.5 

The  data  was  fit  to  an  MNPL  model  that  captured  modal  transition  behavior  and  other 
trends  quite  well,  except  possibly  for  the  out-of-phase  data.  Observed  crack  turning 
angles  show  excellent  correlation  with  the  model,  as  presented  in  Figure  8.  Hollow 
symbols  represent  model  predictions,  and  (neighboring)  solid  symbols  are  test  data,  with 
both  primary  and  secondary  angles  plotted.  Note  that  Highsmith  defined  the  experiments 
with  shear  of  opposite  sign  to  the  data  in  Figure  3,  resulting  in  “flipped”  plots. 
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Figure  2.2  Correlation  of  Mixed-Mode,  Non-Proportionally  Loaded  Specimens  Data  with  Predicted 

Kink  Angles 

The  ability  to  correlate  the  data  was  encouraging,  and  fracture  surface  examination 
(Figures  2.2a  and  2.2b)  showed  a  striking  difference  in  appearance  for  tensile  and  shear 
dominated  fracture,  lending  additional  credibility  to  the  existence  of  a  fracture  mode 
transition.  On  the  down  side,  it  was  found  that  the  data  could  be  fit  fairly  comparably 
using  more  than  one  combination  of  parameters,  which  would  predict  significantly 
different  growth  rates.  Unfortunately,  crack  growth  rates  were  not  measured  during  these 
tests,  so  a  unique  best  fit  could  not  be  established  for  lifing  purposes. 
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a.  In-phase 


b.  Constant  tension/cyclic  torsion 


Figure  2.3  Crack  Face  Appearance  for  Tensile  and  Shear  Dominated  Crack  Growth  [28] 

Discussion 


An  implementation  of  an  MNPL  crack  growth  approach  has  been  presented  capable  of 
analyzing  a  constant  amplitude  linear  cycle  between  two  points.  While  this  is  an 
encouraging  advance,  significant  obstacles  remain.  Vibratory  modes  can  exhibit 
subcycles  involving  all  three  fracture  modes,  particularly  when  neighboring  vibratory 
modes  are  active  simultaneously.  Inevitably,  a  real-world  application  of  the  method  will 
require  a  fully  3D  pairing  method,  which  remains  to  be  defined. 
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Also,  evidence  from  other  authors  [21,27]  suggests  that  addition  of  shear  modes  into  an 
otherwise  mode  I  dominated  cycle  can  in  fact  slow  or  arrest  the  crack  in  some  regimes,  as 
opposed  to  speeding  it  up  as  would  be  predicted  with  positive  fracture  mode  asymmetry 
ratios.  Further,  Tschegg  [20]  showed  that  cracks  propagating  (planar)  in  pure  mode  III 
can  slow  down  and  arrest  even  when  the  applied  AKm  increases  with  crack  length.  As 
aluded  to  earlier,  these  shortcomings  are  likely  attributable  to  roughness  induced  closure 
associated  with  shear  mode  growth,  which  is  not  correctly  accounted  for  by  the  simple 
mode  I  plasticity  induced  closure/friction  method  presented  herein.  It  is  nevertheless 
hoped  that  the  framework  provided  will  remain  useful  as  further  advances  remedy  these 
deficiencies. 
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